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ULTRA-HIGH-ENERGY COSMIC RAY ACCELERATION BY MAGNETIC 
RECONNECTION IN NEWBORN PULSARS 

Elisabete M. de Gouveia Dal Pino ,2 & Alex Lazarian 3 
RESUMEN 

El resumen sera traducido al espanol por los editores. 
We here investigate the possibility that the ultra-high energy cosmic ray (UHECR) 
C> events observed above the GZK limit are mostly protons accelerated in reconnec- 

' tion sites just above the magnetosphere of newborn millisecond pulsars that are 

, originated by accretion induced collapse (AIC). 

X> • ABSTRACT 

["T i ' We here investigate the possibility that the ultra-high energy cosmic ray 

(UHECR) events observed above the GZK limit are mostly protons accelerated 
in reconnection sites just above the magnetosphere of newborn millisecond pulsars 
that are originated by accretion induced collapse (AIC). 
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O !' The detection of cosmic ray events with energies beyond 10 eV (UHECRs) poses a challenge for the 
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understanding of their nature and sources. If UHECRs are mostly protons, then they should be affected by 
the expected Greisen-Zatsepin-Kuzmin (GZK) energy cutoff (~ 5 x 10 19 eV), which is due to photomeson 
production by interactions with the cosmic microwave background radiation, unless they are originated at 
distances closer than about 50 Mpc (e.g., Medina Tanco, de Gouveia Dal Pino & Horvath 1997). On the other 
hand, if the UHECRs are mostly protons from nearby sources (located within ~ 50 Mpc), then the arrival 
directions of the events should point toward their sources since they are expected to be little deflected by the 
intergalactic and Galactic magnetic fields (e.g., Medina Tanco, de Gouveia Dal Pino & Horvath 1998). The 
present data shows no significant large-scale anisotropy in the distribution related to the Galactic disk or the 
local distribution of galaxies, although some clusters of events seem to point to the supergalactic plane (Takeda 
et al. 1999). 

We here discuss a model in which UHECRs are mostly protons accelerated in magnetic reconnection sites 
outside the magnetosphere of very young millisecond pulsars being produced by accretion induced collapse 
(AIC) of a white dwarf (de Gouveia Dal Pino & Lazarian 2000; hereafter GL2000). When a white dwarf 
reaches the critical Chandrasekhar mass ~ 1.4 M Q through mass accretion, in some cases it collapses directly 
to a neutron star instead of exploding into a supernova. The accretion flow spins up the star and confines the 
magnetosphere to a radius Rx where plasma stress in the accretion disk and magnetic stress balance (Arons 
1993). At this radius the equatorial flow will divert into a funnel inflow along the closed field-lines toward the 
star, and a centrifugally driven wind outflow (see Fig. 1 of GL2000). To mediate the field lines of the star with 
those opened by the wind and those trapped by the funnel inflow emanating from the Rx region a surface of 
null poloidal field forms (e.g., Shu et al. 1994). This reconnection region dominated "helmet streamer", will 
release magnetic energy that will accelerate particles to the UHEs. 
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A primary condition on the reconnection region for it to be able to accelerate particles of charge Ze to 
energies E is that its width ARx > 2tl, where Tl is the particle Larmour radius tl = E/ZeBx , and Bx is 
the magnetic field (normal to particle velocity) at the Rx region. This condition and the field geometry imply 
(GL2000): 

b 13 ^ z-i E 2Q n-^ {^f^) 1 ' 2 (i) 

where B i3 is the stellar magnetic field in units of 10 13 G, Q.2.hk is the stellar angular speed in units of 2500 s -1 , 
and E20 is the particle energy in units of 10 20 eV. We find that stellar magnetic fields 10 12 G < B± £ 10 15 G 
and angular speeds 4 x 10 3 s _1 ^ > 10 2 s _1 (or spin periods 1 ms <; < 60 ms), are able to accelerate 
particles to energies P20 J> 1. 

A newborn millisecond pulsar spins down due to magnetic dipole radiation in a time scale given by t* = 
f^/O* ~ 4.3 x 10 7 s B^ n~ 2 5k . We can show that the condition that the magnctosphcre and disk stresses 
are in equilibrium at the inner disk edge results a disk mass accretion rate that is super-Eddington. This 
supercritical accretion will last for ~ r*. As it approaches the end, the newborn pulsar decreases its rotation 
speed due to electromagnetic radiation at a rate r" 1 . The spectrum evolution of the accelerated UHECRs 
is thus determined by r" 1 . The particle spectrum N(E) is obtained from N = N(E) ^§ = N(E) Jf- f2* 
(GL2000): 

N(E) ~ 5.8 x 10 34 GeV- 1 £ Z~^ 2 B u 1/2 E^' 2 ( ^f^ ) ^ (2) 

where £ is the reconnection efficiency factor; the derived spectrum above is very flat which is in agreement with 
the observations. 

The total number of objects formed via AICs in our Galaxy is limited by nucleosynthesis constraints to a 
very small rate ~ 10~ 5 yr^ 1 . Hence, the probability of having UHECR events produced in the Galaxy will be 
only P ~ fbTAic 1 1 — 2 x 10~ 6 , where fb ~ (ARx / Rx) 2 — 10 -2 is the emission beaming factor caused by 
the magnetic field geometry, and t = 20 yr accounts for the time the UHECR events have been collected in 
Earth detectors so far. Since the individual contribution to the observed UHECRs due to AICs in our Galaxy 
is so small we must evaluate the integrated contribution due to AICs from all the galaxies located within a 
volume which is not affected by the GZK effect, i.e., within a radius R^q — i?c/50 Mpc. Assuming that 
each galaxy has essentially the same rate of AICs as our Galaxy and taking the standard galaxy distribution 
n G ~ 0.01e ±a4 /i 3 Mpc -3 (with the Hubble parameter defined as H a = h 100 km s _1 Mpc -1 ), the resulting 
flux at E 2 q > 1 is F(E) ~ N(E) n G t A i C ~ y R g , which gives 

F{E) ~ 1.1 x 10- 2 ^GeV- 1 cm- 2 S - 1 Z-^B^ /2 E w 3/2 r AIC , 5 - 1 n OM R 50 ( ^q^ )"^ (3) 

where taic,^ 1 = taic^ 1 /10~ 5 yr -1 , and no. 01 = nc/0.01 h 3 Mpc -3 . Observed data by the AGASA exper- 
iment (Takcda et al. 1999) gives a flux at E =10 20 cV of F(E) ~ 4 x 10 -30 Gcv -1 cm -2 s" 1 , so that the 
reconnection efficiency factor needs to be only £ ^ 3.6 x 10~ 3 in order to reproduce the observed signal. 
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